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1 Introduction 



Lattice gauge theories ||T| are the only known rigorous non-perturbative formulations of 
non Abelian theories and their interpretation as regularized versions of continuum quantum 
field theories in Euclidean space has been among the most fruitful theoretical ideas in the 
last years. They also offer the unique possibility to compute correlation functions non- 
perturbatively through numerical simulations, and therefore they represent a formidable 
laboratory where formal propositions may be tested and fundamental phenomenological 
quantities can be computed from first principles. 

Lattice gauge theories are defined on a discretized space-time which cuts-off the high and 
low frequences and renders the theory finite. The fundamental gauge fields are elements of 
the underlying group and, since the group is compact, the lattice functional integrals are 
well defined without any gauge-fixing. 

In the following we mostly consider pure Yang-Mills theories (without fermions). Yet it is 
interesting to note that, if the fermionic operator satisfies the Ginsparg- Wilson relation 0, 
fermions can be introduced on the lattice preserving chiral and flavor symmetries at finite 
cutoff. 

Non-perturbative lattice gauge- fixing becomes unavoidable to extract information from gauge- 
dependent correlators . It is necessary in order to study the propagators of the fundamental 
fields appearing in the QCD Lagrangian in the non-perturbative region. 

It is also necessary in some non-perturbative renormalization schemes @, ^ which use gauge 
dependent matrix elements to renormalize composite operators, and it can become a funda- 
mental technical ingredient in the so called non gauge invariant quantizations of chiral gauge 
theories These motivations justify the efforts to obtain a consistent non-perturbative 
lattice gauge fixing. 

In this review we will discuss some of the problems in lattice gauge fixing, selecting the topics 
on the basis of their importance and of our personal experience in the field. 

In Sec. ^ we review the most popular approach to define gauge-fixing on the lattice: we give 
the definitions of the gauge dependent correlation functions and we sketch briefly the steps 
of the numerical procedures adopted. 

In Sec. 1^ we review the non-perturbative definition of the Landau and Coulomb gauges and 
we describe the main algorithms used in the literature to enforce these gauges numerically. 

In Sees ^, ^, 1^, we review other lattice gauge conditions and the corresponding gauge- fixing 
procedures proposed in the literature. 

In Sec. ^ we briefly sketch the role of the gauge choice in understanding the physics of quark 
confinement. 

Sec. H is devoted to the gauge fixing implementation in the Langevin dynamics algorithm, 
which is necessary to overcome divergent fluctuations along the gauge directions. 

In Sec. |10] we discuss some problems related to the ambiguities in the lattice definition of 
the gauge potential. Different regularized definitions and their effects on gauge dependent 
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quantities are analyzed. 

In Sec. ^ we discuss the problem of numerical Gribov copies, their effects on physical 
quantities and some approaches which have been proposed to remove this ambiguity. 



Sec. |T2| is devoted to the lattice QCD gauge dependent smoothing procedures. 



In Sec. |I3] we draw our conclusions and acknowledgements. 



2 Standard Non-Perturbative Gauge Fixing 

In the standard formulation of lattice gauge theories proposed by Wilson [0 the link f/^(x) 
are the fundamental gauge fields of the theory, they are group elements of SU (N) in the fun- 
damental (N-dimensional) representation and they transform under a gauge transformation 

G{x) as 

U^ix) = Gix)U^ix)G\x + fi). (1) 
The gauge invariant action is defined as 

-5 = [l - ^Tr [P,^{x) + Plix)] ] (2) 

plaq 

where, in the standard notation, /3 = 2N/gl, go is the bare coupling constant, Pfj,u{x) is the 
Wilson plaquette, i.e. the path-ordered product of link variables 

P,u{x) = U,{x)U,{x + fi)Ul{x + u)Ul{x) (3) 

around the boundary of a plaquette P. The expectation value of any gauge invariant operator 
0{U) is given by 

{O) = ^ J dUe~^^^^0{U) (4) 



where dU denotes the group-invariant integration measure over the links satisfying the fol- 
lowing properties: 

j h{U)dU = j h{VU)dU = j h{VU)dU V 1/, f/ G G (5) 

dU = I- 

being h{U) a generic function of the links. Since the domain of the link integration is com- 
pact, the lattice functional integrals (H) are well defined and the gauge invariant correlation 
functions can be computed without fixing the gauge. 
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In the ideal case where the gauge-fixing condition f{U^) = has an unique solution for 
each gauge orbit, i.e. there are no Gribov copies [0, the Faddeev-Popov procedure |^ can 
be applied. The gauge-invariant Faddeev-Popov determinant A/(f/) is defined as 

Af{U) I dG5U\U'')) = 1 , (6) 

where the integration is over all gauge transformation G. By inserting the previous identity 
in the functional integrals (^), changing the variables U —>■ U'^ and the order of integration, 
using the gauge invariance of dU, S{U), Af and of 0{U), we can write 

{0)f = \j dUe-'^^^Af{U)6{f{U))0{U) . (7) 

For gauge-invariant quantities this expression is equivalent to Eq. (H). Eq. (|^) is the Faddeev- 
Popov definition of the correlation functions of gauge- dependent operators. 

In absence of Gribov copies, the Faddeev-Popov determinant can be expressed as an integral 
over the ghosts and anti-ghost fields 1] and f/, obtaining 

_ JdUd\df]dri e-^(^)e-ti'^'e'^J^^C'([/) 
^ fdUd\dr]dr]e-^(u^e-^^^'e^ff^ 

where A are Lagrangian multipliers and 6 represents the lattice BRST transformations 
defined as 

5f] = iX (9) 

6X = . (10) 

The gauge-fixed action in (|^), including the ghost terms, is local and is invariant under 
BRST transformations. The Faddeev-Popov procedure can be replaced by the more formal 
apparatus of the BRST symmetry. It resembles the continuum formula 



(C») = j M^5r75r/Ce-^(^)~^«'^-'('''*''^)5(/(A)) . , (11) 

On the contrary in presence of Gribov copies, i.e. multiple solutions of the equation 

/([/«) = (12) 

for a given gauge configuration U, the Faddeev-Popov determinant cannot be expressed as 
an integral over the ghosts and the BRST invariance is lost. Labelling the different solutions 
of (|T^ by Gi, in the integral functional (|^) we sum over several gauge equivalent copies of 
the same configuration obtaining 



Idet |g=g. 
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Eq. (0) leads to an acceptable but very inconvenient gauge fixing procedure. 

Alternative procedure to mantain the BRST symmetry in the gauge fixing process despite 



of the presence of the Gribov copies have been proposed . But Neuberger has shown |[TT[] 
that on the lattice the requirement of the standard BRST invariance of the gauge-fixed 
action leads to the non-perturbative level to disastrous results, i.e. the physical observables 
are reduced to the an undetermined form. The argument can be sketched as follows : let 
us define the function 



Fo{t) 



dU dXdrjdr] e 



which satisfy 



dFo{t) 
dt 



dUdXdrjdf] 6 







(14) 



(15) 



because the integral of a total BRST variation vanishes identically. On the other hand 
-^0(0) = because the integrand does not contain ghosts and therefore Fo{l) = and 
{O) = ^. Possible ways out of this paradox have been proposed in O, O by imposing a 



y Possible ways out of this paradox have been proposed in 
modified BRST symmetry on the lattice which converges towards the conventional one in 
the continuum limit, interesting remarks on this subject can also be found in [Q. 



The standard numerical gauge-fixing procedure on the lattice |]T5[ is obtained by reversing 
the Faddeev-Popov analysis described above and neglecting the presence of Gribov copies. 
In (^, by multipling by J dG (nothing depends on G), changing the order of integration, 
changing variables U W^, using the gauge invariance (^ of dU, S{U), and performing 
the integral over G (taking into account the 6 function) we can write 



(16) 



where G{U) is the gauge transformation for which f{U'^^'^^) = 0. The definitions in Eqs. (|^ 
and (|1^) of the correlation functions of gauge dependent operators are equivalent if the 
equation fiU'^) = has a unique solution, i.e. there are no Gribov copies. In this case 
the Faddeev-Popov determinant A/(f/) cancels out when the integral of the 6 function is 
performed. 

The Eq. (^) summarizes the procedure implemented by the numerical algorithm |T^: 



a set of thermalized configurations {C^} is generated with periodic boundary con- 
ditions according to the gauge invariant weight e""^*^^^; 

for each {Cu} a numerical algorithm computes the gauge transformation {G}; 

each thermalized configuration {Cu} is gauge rotated obtaining the set of the gauge 
fixed thermalized configurations of the gauge-fixed set {C^g}; 

the expectation value of an operator is given by the average of the values of the operator 
evaluated at each gauge rotated configuration of the gauge-fixed set: 



5 



(O) 



Latt 



1 

N 



{con/} 

The complexity of the ghost technique is replaced by the large amount of computer time 
spent to obtain numerically the gauge transformations which satisfy the gauge condition 
required. The entire procedure is rigorous only when the gauge fixing condition is free 
from Gribov copies. Otherwise the effect of Gribov copies must be taken into account (see 
Section |TT|), since the definition of the correlation functions would depend on the way the 
gauge fixing algorithm selects a preferred Gribov copy. 

In the case of an imperfect or inadequate gauge fixing the measurement of a gauge dependent 
operator is at best affected by additional fluctuations to be summed up to the intrinsic 

In other cases, as in some calculations on U{1) | 



statistical noise 16 



18, 19] and on 



confinement vortex picture [20 



the influence of lattice Gribov copies can mask the 



regular behaviour of a measurement (see later) . 



3 Landau and Coulomb Lattice Gauge-Fixing 



The standard way of fixing the Coulomb and Landau gauges on the lattice |23|, |T^, |2^, [25 
is based on the numerical minimization of the functional 



Fu[G] = -ReTr t/^.^^^^ 



X) 



(18) 



X /i = l 



where / is 3 for Coulomb and 4 for Landau gauge. -^[/[G'] is constructed in such a way that 
its extrema G* 

=0 (19) 

G* 



SG 



are the gauge fixing transformations corresponding to the discretized gauge condition 

A^(x) = J2 (^?(^) - ^?(^ - A)) = , (20) 



where 



U,{x)~Ul{x) 
2iago 



(21) 

Traceless 

The "standard definition" (^) of the gauge potential is naively suggested by the interpreta- 
tion of Ufj_{x) as the lattice parallel transport operator and by its formal expression in terms 
of the "continuum" gauge field variables, 



U^{x) = e 



iagoAfiix) 



:22) 



where a is the lattice spacing. The second variation of Fcoui is the lattice Faddeev-Popov 
operator for the Landau and Coulomb gauge. 
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The lattice gauge-fixing sketched above is the analogous of the continuum Gribov's proce- 
dure. Moreover there is also a corrispondence between the lattice functional in Eq.(]TE|) and 
the continuum one 



iGl |2 



(23) 



Fa[G] = -Tr j d^x (Aj(a;)Aj(a;)) = - {A^,A^) = 

which reaches its extrema when the gauge-fixing condition 

dfj_A^ = 0, supplemented with periodic boundary condition; (24) 

is satisfied, fi goes from 1 to 3 or 4 in the case of Coulomb and Landau gauge respectively. 
It is remarkable that Eq. ([18|) does not correspond to the natural discretization of the con- 



tinuum functional ( p3D according to the lattice definition of the gluon field (^) but it differs 
from that by 0{a) terms. The form in Eq. ([T8|) is adopted not only for its simplicity but 
also because it leads to the gauge condition (PP]). 

The most naive algorithm to minimize the functional in Eq. (^) sweeps the lattice by 
imposing the minimization requirement one site at a time, and repeating the process until 



the gauge transform has relaxated sufficiently into a minimum p2| , [15] . Actually in order to 
fix the gauge one needs just to reach any stationary point of F, hence the requirement of a 
minimum is a somewhat stronger request naturally adopted by the numerical procedure. 

In order to study the approach of the functional ([18|) to a minimum, the values of two 
quantities are usually numerically monitored. The first one is -F1/[G'] itself, which decreases 
monotonically and eventually reaches a plateau. The other one, denoted by 6, is defined as 
follows: 

^""^V = 7 [A^(a;)(A^)t(a;)],:^ J d'x Tr{d,A^^)\ (25) 

X X 

where V is the lattice volume. The function ^ is a measure of the first derivative of -F[/[G'] 
during the gauge-fixing process; it decreases (not strictly monotonically) approaching zero 
when F[?7*^] reaches its minimum. The desired gauge fixing quality is determined by stopping 
the computer code when 9'^ has achieved a preassigned value close to zero which is often 
defined the gauge fixing quality factor. In Fig. |l] it is shown the typical behaviour of 6 and 
of the quantity = |F — Fmin\ as function of the gauge fixing sweeps. These behaviours 
can change quite a lot among different thermalized configurations. 

The choice of the gauge fixing quality is a delicate point in the case of a simulation with 
a large volume and a high number of thermalized configurations. Of course, the better is 
the gauge fixing quality, the more computer time is needed. Moreover it is impossible to 
know in advance, before computing the gauge dependent correlation functions, whether an 
a priori criterion is suitable or not. So that, the stopping 9 value is normally fixed on the 
basis of a compromise between the available computer time and the gauge fixing quality. 
Sometimes, in the case of calculations performed on computers with single precision fioating 
point, the maximum gauge fixing quality is limited by a value of the order of the fioating point 
zero: 6 ~ 10~^, this value is usually enough to guarantee the stability of gauge dependent 
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Figure 1: Typical behaviour of 6 (curve a) and Fr (curve b) as function of the number of 
gauge fixing sweeps for the Landau gauge. Lattice size is 8^ • 16. 
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correlators even in high precision measurements, hke the calculations of gluon and quark 
propagators (see for example Ref. Q]) and the calculus of the running QCD coupling by 
means of the 3-gluon vertex function. In these cases gauge fixing becomes a time consuming 
part of the computation, comparable to the calculation of a quark propagator. 

The lattice Landau and Coulomb gauge-fixing are affected by the problem of lattice Gribov 
copies. Section ^ is devoted to this issue. 

It is also interesting to see the gauge-fixing procedure from another point of view: it can 
be considered as the process of finding the ground state of a dynamical system (a spin 
system like in the Ising model) where F takes the place of the hamiltonian, the gauge 
transformations G"s are dynamical variables, belonging to the SU{3) group, and the links 
are the couplings. This analogy is clearly seen by writing the gauge trasformation (|T|) in the 
functional form (PT 



3.1 Acceleration of Landau and Coulomb Gauge-Fixing 



For large lattices the gauge fixing algorithms, as other iterative methods, converge slowly 
due to long-range correlations and large condition numbers of the matrices which control 
the algorithms. This is a crucial problem usually called critical slowing- down (see for ex- 
ample [^). To reduce the critical slowing-down in gauge fixing algorithms, two classes of 
improvements are often adopted: 



overrelaxation, originally proposed in Ref. |28|, to speed up gauge-fixing algorithms 



Fourier preconditioning to adjust the matrix governing the system evolution so that 
all the eigenvalues become approximately equal to the largest one without affecting the 
final answer. 



The overrelaxation algorithm, a technique originally introduced to improve the convergence 
of iterative methods to solve classical linear algebra and differential equation problems, is 
particularly suitable to face critical slowing down in numerical simulations, as first shown 
in Ref. |2^. The effectiveness of this method has been studied in different papers; for a 
general discussion see also Refs. |5D|, ^ |53|, 0, ^ |5B[ . The overrelaxation method is 
implemented in the process of gauge fixing by replacing G{x) with its power G'^(a;), at each 
iteration. In practical computations, C^lx) is given by a truncated binomial expansion 



N 

E 

n=0 



{G-I)\ 7n(^) 



r(c^ + ii 



r(cu + i 



n] 



(26) 



where 1 < u < 2 and usually 2 < < 4. The u parameter is tuned empirically at an 
optimal value ujopt, (a typical value is Uopt — 1-75), to reach the fastest convergence. Before 
G"^(x) is applied on the link, it has to be appropriately normalized to belong to the gauge 
group. 

The study of the convergence as a function of the overrelaxation parameter is shown in Fig. |^. 



It shows clearly that it is possible to distinguish two different regimes [28, 29, 30, 36]. The 
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Figure 2: Efficacy of gauge fixing as a function on the overraxation parameter uj from 
Ref. The plot shows the decrease of 9 for different values of uj and a 4^ lattice at 

(5 = 5.6. 



initial behaviour appears to be almost insensitive to the value of u and corresponds to large 
local fluctuations from one iteration to another. In the second regime the fluctuations are 
smooth, characterized by a long-range pattern, and the rate of convergence is quite sensitive 
to the u value since here one faces the problem of critical slowing down and the optimization 
of the algorithm becomes crucial. On the lattices considered in Ref. a factor of 3 or 
5 is gained in the number of iterations while a computer time overhead of a factor 1.5 for 
iteration is paid so that the resulting gain in time due to this procedure is about a factor 
of 2. Another small price to pay is related to the need for tuning empirically u at an optimal 
value LJopt for which convergence is most rapid, fortunately this value does not depend too 
much on different configurations at fixed (3 and V. It is remarkable that any algorithm 
for Landau or Coulomb gauge fixing can be easily modified in such a way to include the 
overrelaxation, the only change required is to add the expansion (^61) at the end of each 
iteration of the normal gauge fixing procedure algorithm. A variation of the overrelaxation 
method is the stochastic overrelaxation algorithm, proposed by Ref. [^, in which a local 
gauge transformation G'(x)^ replace G{x) with probability p. The actual acceleration gain 
turns out to depend strongly on p and the procedure definitely diverges for p = 1. 

Monitoring the deviation from the Landau gauge at each site, it is possible to see a broad 
range of small deviations while in the Fourier space many slowly decaying modes are ob- 
served, including, but not limited to, the longest wavelenghts. 

The Fourier acceleration (FA) technique applied to the gauge fixing, alleviates the problem 
of critical slowing down. The idea is to precondition the problem using a diagonal matrix in 



momentum space which is related to the solution of a simplified version of the problem El, 
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• The overall relaxation time is determined by the smallest eigenvalue in the momentum 
space of the matrix governing the iterative algorithm. 

Here we describe the use of Fourier acceleration to improve the convergence of the the Landau 



gauge fixing algorithm, following Ref. pj]. This method modifies the gauge transformation 



G{x) in Fourier space in such a way that all modes converge as fast as the fastest mode. 
In order to fix a lattice version of Landau gauge one minimizes the expression (|18D in the 
space of gauge-equivalent field f/^*^ defined in relation (|I]). Following a naive steepest-descent 
method, we differentiate with respect to the gauge trasformation, and at each step of the 
iterative procedure, G{x) takes the expression: 

G{x) = exp{^[j2 [A-.(f/.(x) - Ulix)) - l/iVcTr[A_,(f/,(x) - Ul{x))]\]] (27) 

V 

where A,^(f/^(x)) = U^{x — v) — U^{x) and a is a tuning parameter. To optimize the 
convergence, the Fourier accelerated method replaces the Eq. (E^) by 



2 2 

G{x) = exp[F-^(|^i^f^F(5^ A_,(f/,(x) - U\x)) ~ trace^^] (2 



,2 

where are the eigenvalues of the lattice version of the (9^) operator, a is the lattice spacing 
and F denotes the Fourier transform operator. Thus, in this case, the preconditioning is 
obtained using in momemtum space a diagonal matrix with elements given by . 



The FA method in its original form makes use of the fast Fourier transform algorithm to 
evaluate F and which requires a computer time proportional to VlogV where V is 

the lattice volume p4|, making it very appealing from the numerical point of view. Note, 
however, that on parallel machines the cost of the Fourier transform is not negligible because 
of its high non-locality. For this reason the FA is not considered appealing anymore in 
large scale simulations. An implementation of the FA method for Landau gauge fixing, 
avoiding completely the use of the Fourier transform, has been proposed and tested for the 
4-dimensional SU{2) case, on serial and on parallel machines, in Refs. ||38|, In Figs. 



and ^ the behaviours of a gauge fixing evolution with and without the Fourier accelerated 
algorithm as a function of the iteration number are shown. 



In Ref. global gauge fixing on the lattice specifically to the Landau gauge, is discussed 
with the goal of understanding the question of why the process becomes extremely slow 
for large lattices. The author constructs an artificial "gauge fixing" problem which has 
the essential features encountered in the real case. In the limit in which the size of the 
system to be gauge fixed becomes infinite, the problem becomes equivalent to finding a series 
expansion in functions which are related to the Jacobi polynomials. The series converges 
slowly, as expected. It also converges non-uniformly, which is an observed characteristic 
of gauge fixing. In the limiting example of Ref. |^ the non-uniformity arises through the 
Gibbs phenomenon. 

Gauge fixing algorithms have also been used as a prerequisite for the Fourier acceleration of 
other procedures required for simulations of lattice gauge theories as the matrix inversion, 
for algorithms based on the Langevin equation and others. 
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Figure 3: 9 plotted as a function of iteration number for gauge fixing with and without 
Fourier acceleration, redrawn from Ref. |23 . 




200 



Figure 4: Fjj (only the first 200 steps are shown) plotted as a function of the gauge 
fixing iteration number for the same configuration as in Fig. ^ with and without Fourier 



acceleration, redrawn from Ref. p3 
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A review of the performances of some standard gauge procedures for non-Abelian gauge 
theories can be found in Ref. 



Techniques involved in gauge fixing, as reunitarization 
and convergence criteria, accelerating procedures and performance of algorithms on the 
parallel machines CM2 and CMS are also analyzed. Critical slowing-down of several gauge 
fixing algorithhms for different gauges in the SU{2) case, at zero temperature, is discussed 
in Ref. fHl 



4 The Soft Covariant Gauge 

A nonperturbative method for gauge fixing in the continuum has been proposed by Jona- 
Lasinio and Parrinello and by Zwanziger and extended on the lattice in |Q. They 



suggest to modify the gauge invariant Wilson's partition function 

Z = j dU e'^^^\ (29) 
by simply inserting an identity in it (^9]) : 

Z^,, ^ jdU e-^(^) r'[U] J dG e-/5^^'^I^''l. (30) 

F\U'^] is a generic function of the links not invariant under general gauge transformations 
and the gauge invariant quantity I[U] is given by: 

I[U] = [dG e-/5*^'^[^^l. (31) 



Eq. (^OD corresponds to the first step of the standard Faddeev-Popov gauge-fixing procedure. 
However, unlike what would happen in the continuum, Zmod = Z is a finite quantity because 
the group of gauge transformations on a finite lattice is compact. Thus Zmod can provide a 
new definition for the expectation value of 0[f/]: 

< O >w= / dU e-^(^) r'[U] I dG e-^^^^^l^"] 0[U% (32) 

If 0[U] is a gauge invariant operator then < O >=< O >mod, while eq. (^) defines the 
expectation value of gauge-dependent operators. 

By defining 

< 0[U] >G= I-'[U] j dG e-'^M^^if/G] Q^^G^^ ^33) 
< O >mod can be cast in the form: 



JdU e-^(^) <0[U]> 



< O >raod= ' \ ^JnC'^'' =< <0[U]>G>. (34) 



JdU 



-S(C/) 
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The above expression indicates that in the gauge-fixed model the expectation value of a 
gauge dependent quantity 0\U] is obtained in two steps. First one associates with 0\U] the 
gauge invariant function < 0\U] >g, which has the form of a Gibbs average of ©[t/*^] over 
the group of gauge transformations, with a statistical weig ht factor e-^^'^^^ I Then one 
takes the average of < 0[U] >g a la Wilson. 

This suggests the following numerical algorithm: 

• generate a set of link configurations Ui, . . . Un, weighted by the Wilson action, via the 
usual gauge invariant Monte Carlo algorithm for some value of f3; 

• use each of the Ui as a set of quenched bonds in a new Monte Carlo process, where 
the dynamical variables are the local gauge group elements G{x) located on the lattice 
sites. These are coupled through the links Ui, according to the effective Hamiltonian 
F[Uf]. In this way one can produce for every link configuration Ui an ensemble of 
gauge-related configurations, weighted by the Boltzmann factor exp{—l3M'^F[U^]). 
We call < 0[Ui] >g the average of a gauge dependent observable O with respect to 
such an ensemble, in the spirit of Eq. (|33). 



Finally, the expectation value < O >mod is simply obtained from the Wilson average 
of the < 0[Ui] >G, i-G.: 



1 ^ 



i=l 



In the above scheme can be interpreted as a gauge parameter, which determines the 
effective temperature 1//3M^ of the Monte Carlo simulation over the group of gauge trans- 
formations. 



Adopting as gauge fixing action the form of the Landau gauge fixing functional (|T8D it is 
possible to get the connection between this scheme and the usual Landau gauge fixing. It 
turns out that the stationary points of F[f/*^] correspond to link configurations U^ that 
satisfy the lattice version of the Landau gauge condition. All such configurations correspond 
to Gribov copies. In particular, those corresponding to local minima of F[f/*^] also satisfy 



a positivity condition for the lattice Faddeev- Popov operator [^. As a consequence, in 
the limit oo, the above gauge- fixing is equivalent to the so-called minimal Landau 

gauge condition, which prescribes to pick up on every gauge orbit the field configuration 
corresponding to the absolute minimum of F[f/*^] ||47|| . 

While this method is conceptually very simple, much less is known about its perturbation 
theory expansion, as compared with the standard gauge fixing based on BRST invariance. 
In fact in the Faddeev-Popov case the determinant can be expressed as an integral over 
ghost fields, and the gauge-fixed action including the ghost terms is local, whereas here 
this is not the case. This is an important difference, because locality is a key ingredient 
in power-counting arguments, and thus at the heart of the usual perturbative analysis of 
renormalization. 
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It is therefore of interest to find out whether perturbation theory can be systematically 
developed for gauge-fixed Yang-Mills theory correlated with the soft covariant gauge fixing, 
and its relation to the usual Faddeev-Popov procedure. This question has been addressed 
some time ago in Ref. and recently in Ref. \W, 0, A numerical implementation 

of this technique has been applied |Q to a study of the gluon propagator on small lattice 
volumes but the applicability of this method to physical lattices seems to be numerically 
demanding. 



5 The Generic Covariant Gauge 



In the continuum 53 1, the Faddeev-Popov quantization for covariant gauges is obtained by 



fixing the gauge condition 



d.A'iix) = A{x) (36) 



where A{x) belongs to the Lie algebra of the group. Since gauge-invariant quantities are not 
sensitive to changes of gauge condition, it is possible to average over A(a;) with a Gaussian 
weight. As usual the Faddeev-Popov factor can be written as a Gaussian integral of local 
Grassman variables, the resulting effective action is invariant under the BRST transforma- 
tions and the correlation functions of the operators satisfy the appropriate Slavnov- Taylor 
identities. The expectation value of an operator can be cast into the following form (to be 
compared with Eq. (p^)): 

(O) = J 5Ae-^^'^*"^"(^') J 5A^5r/5r7Ce-^(^)-^«'^-*(^'^'^)5(9^A^ - A) , (37) 

obtaining 



In the perturbative region, the renormalized correlation functions can be compared with the 
same quantities computed in the standard perturbation theory. 



In Ref. [54 1 it has been proposed a numerical procedure to implement this covariant gauge- 
fixing on the lattice. The algorithm is based, as in the Landau case, on the minimization 
of a functional Ha[G\ chosen in such a way |^6[ that its absolute minima correspond to a 



gauge transformation G satisfying the general covariant gauge-fixing condition in Eq. (pG]). 
In order to fix the covariant gauge on the lattice, one should be able to find a functional 
H{G) stationary when the gauge (|36|) is fixed. The most simple way to define H{G) would 
be to find a functional h{G), to be added to F, in such a way that 

5G - ^ " 



It has been shown [HR that this functional does not exist for a non abelian gauge theory. 



It is interesting to give the outline of the proof. Writing the gauge transformation: 

G = exp(i^w"T") 
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where T'^ are the eight SU{3) generators, the first derivative of F takes this form showing 
the stationarity of F{G) when 9^^^ = 0: 



SF{G) 
5w^ 



where 



e-f - I 

7 



-1 ab 



1 



ab fabc„,,c 



Then the derivative of h should have this form: 



Sh 



- (A"(x)$"^(w(x))) . 



6w^{x) w 

A necessary condition for the existence of such a functional would be 

5^h 5^h 



(39) 



(40) 



(41) 



5w^{x)5w^{y) 5w^{y)5w^{x) 
which implies the integrability condition 



(42) 



6 



6 



(A"(x)<l>""(ti;)) 



(43) 



Expanding $"^(w(x)) in power of w{x), the Eq. (|43|) should be satisfied order by order in 
w{x). From Eq. (BO) one has 



(44) 



Equation ( ^3] ) is then in contrast with the antisymmetry of f"'^'^. 
The new functional, chosen to resemble F for this gauge, is ^6 



Ha[G] = J d^xTi [{d^A<; - A)(9,A^ - A 



(45) 



which obviously reaches its absolute minima {Ha[G] = 0) when Eq. ( p^ ) is satisfied. There- 
fore in this case the Gribov copies of the Eq. (|36|) are associated with different absolute 
minima of Eq. (|^) . Due to the complexity of the functional (|45| ) , it may also have relative 
minima which do not satisfy the gauge condition in eq. (|36D (spurious solutions) |^ because 
the stationary points of Ha[G] actually correspond to the following gauge condition 



(46) 



where D^, is the covariant derivative. The spurious solutions, therefore, correspond to zero 
modes of the operator D^^dy. Of course the numerical minimization of the discretized version 
of eq. ( ^5l) can reach relative minima (spurious solutions) with Ha[G\ which are not dis- 
tinguishable from the absolute minima. Hence this could simulates the effect of an enlarged 
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set of numerical Gribov copies. Preliminary checks at a = do not show any practical 
difference between the use of the new functional with respect to the standard Landau one. 

On the lattice, the expectation value of a gauge dependent operator O in a generic covariant 
gauge is 



-(3S{U) 



(47) 



that is the straightforward discretization of Eq. (|37| ) where Ga is the gauge transformation 



that minimizes the discretized version of the functional (|45|) . In order to avoid a quadratic 
dependence on G of Ha[G] during the single, local minimization step of the gauge fixing 
algorithm, the discretization of Ha[G] has been done by modifying, in each different term of 
Hu[G], the definition of A by terms of order a, "driven discretization". The proposed form 
of Hu[G] on the lattice is the following: 



Hu[G] 



VTa^g 



;Tr^j'^{x)j'^\x) 



(4J 



where 



J{x) = N{x) - igA{x) , 
N{x) = -8I + J2{UKx-u) + U,{x)) 



(49) 



Hu[G] is positive semidefinite and, unlike the Landau case, it is not invariant under global 
gauge transformations. The functional Hu[G] can be minimized using the same numerical 
technique adopted in the Landau case. In order to study the convergence of the algorithm, 
two quantities can be monitored as a function of the number of iteration steps: Hu[G] itself 
and 

Oh = T^J^^rlA^K] , (50) 



where 



Affix) 



Xh{x) 



oc 



SHu[G] 



Traceless 



6e 



(51) 



and 



Xnix) 



= ^(f/^(a;)J(a; + /x) + f/i(x-/i)J(x-/i)) 
- 8J{x) -721 + igN{x)A{x) . 



(52) 



Ah is the driven discretization of the eq. (|46D supplemented with periodic boundary condi- 
tions; it is proportional to the first derivative of Hu[G] and, analogously to the continuum, 
it is invariant under the transformations A{x) A{x) + G, where C is a constant matrix 
belonging to the SU{3) algebra. During the minimization process Oh decreases to zero and 
Hu[G] becomes constant. The quality of the convergence is measured by the final value of 
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Figure 5: Comparison of the behaviour of the gluon propagator transverse part as function of 
lattice time, at two different values of the gauge fixing parameter: a = (corresponding to the 
Landau gauge), and a = 8, for a set of 221, SU{3) configurations at /? = 6 with volume=16^ • 32. 
Prom Ref. |5^. 



9h- We refer to the original paper \M\ for the discussion of consistency checks and further 



details. In Fig. |^ we plot the behaviours of the gluon propagator tranverse part measured 



using this technique to fix the gauge at two different values of the gauge parameter a ^1 . 
A sensitive dependence of the gluon propagator transverse part on the gauge parameter is 
clearly reported. The simulation has been performed over an ensemble of 221, j3 = 6 SU{3) 
thermalized configurations with volume 16^ • 32. and has turned out to be moderately time 
consuming. 

Many interesting considerations on the gauge fixing related to the gluon propagator can 
be found in the review 0, and about the relationship between the gluon propagator and 



confinement, in two recent papers |59[]. 



6 The Laplacian Gauge 

The Laplacian gauge was proposed in alternative to the standard gauge-fixing procedures 
in order to have a smooth gauge- fixing which overcome the problem of Gribov's ambiguities 
^|. The smooth configuration is obtained by rotating the gauge in such a way that the 
eigenvectors corresponding to the smallest eigenvalues of the covariant Laplacian are smooth 
functions of the lattice coordinates. The lattice covariant Laplacian is defined as 

A{U)ab{x, y) := ^[25{x - y)5ab - U^{x)abS{x + fi~y)- U^{x)lf^6{y + p, - x)] , (53) 
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and its eigenfunctions are defined by: 



AiUU{x,y)fay) = >^'f:{x) (54) 

where > are the eigenvalues. We have suppressed the gauge field indices on U{x) G 
SU{N). The gauge transformation G{x) that defines the Laplacian gauge is computed from 
the eigenfunctions with the lowest eigenvalues in order to select the smooth modes in 
the gauge field. 

Specializing to gauge group SU{2) the eigenvalues have a twofold degeneracy, due to the 
charge conjugation symmetry U = a2U*a2, f'^ — > a2f^* ■ The au are the usual Pauli matrices. 
The two degenerate eigenfunctions with the smallest eigenvalue, and cr2/°*, define a 2 x 2 
system on all sites x, which is projected on SU{2) to obtain the gauge tranformation G(x), 

where p(x) = (|/i (a;)P + I/2 (3^)1^)"'^^^ = 1 ^"^^ the two degenerate eigenfunctions /° and o"2/°* 
are normalized, + l/2(^)P) — 1 orthogonal. 

A detailed discussion on the G(c) definition ambiguities rising from lowest eigenvalues degen- 
eration and p(a;)'s zero values can be found in [^]. Nevertheless their effects can be controlled 



by increasing the numerical precision with which the lowest eigenvalues and eigenfunctions 
of the Laplacian are computed. 



The Laplacian gauge on the lattice is investigated numerically in Ref. |6T[] using the gauge 
fields U{1) in two dimensions and SU{2) in four dimensions. The Gribov problem is ad- 
dressed and to asses the smoothness of the gauge field configurations they are compared to 
configurations fixed to the Landau gauge. This comparison indicates that Laplacian gauge 
fixing works well in practice and can offer a viable alternative to Landau gauge fixing. The 
implementation of this gauge for the SU{3) group has been studied in Ref. [|62| . 



A perturbative formulation of the Laplacian gauge for the SU{2) group is presented in 
Ref. however the renormalizability is still to be demonstrated. 



7 The Quasi- Temporal Gauge 



In many cases the Landau or Coulomb gauge fixing consume a large fraction of the compu- 
tational cost of a simulation. Therefore it could be extremely advantageous to find low-cost 
alternative gauges with the features of smoothness and limit to the continuum required by 



the simulations. To this aim, in Ref. 0] a lattice version of the quasi-temporal gauge (QT 
gauge), proposed in Ref. 



35| and formulated rigorously in Ref. has been studied. This 
gauge is a variation of the temporal gauge, widely studied in Ref. ||6^, |6^, |6^. It is defined 
by enforcing the Coulomb condition at a given time t = to: 



d- A {x,to) = V(f,to)- 



(56) 
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together with the temporal gauge: 

= Vx, (57) 

at all points. The association of the Coulomb gauge on one time slice and the axial gauge 
makes the quasi-temporal gauge a complete gauge with the same properties of the Coulomb 
one but with the advantage of being roughly T times cheaper to implement on the lattice 
(T being the time direction length of the lattice). The well-known pathologies of the pure 
temporal gauge are overcome by the trick of time slice fixed into the Coulomb gauge. In par- 
ticular the Gauss' law is satisfied and the problematic pole of the tree level gluon propagator 
is removed. The algorithmical implementation of the quasi-temporal gauge is straightfor- 
ward. Once the Coulomb condition holds at t = to, the temporal gauge (|57|) can be trivially 
imposed by visiting sequentially each timeslice and gauge-transforming the temporal links 
f/o(x, t) into the unit group element. On a periodic lattice, this can be done for all but one 
time tf, so that the temporal links ?7o(x, tj), rather than being unity, end up carrying the 
value of the Polyakov loop at x. Since the computational cost of temporal gauge fixing is 
negligible, it follows that the QT gauge is roughly T times faster to implement than the 
Coulomb gauge. A possible drawback of this gauge condition is that it is not invariant under 
time translation, because of the Coulomb condition at to- Moreover, as well as the Coulomb 
condition, it is affected by the Gribov ambiguity. 

In order to test the feasibility of lattice non-perturbative calculations in the QT gauge, it has 
been calculated the renormalisation constant Za of the axial current by using WI's on quark 
states. This quantity was already measured with several different methods and therefore it 
is useful to test the quasi-temporal gauge. The final numerical results of a simulation on a 
volume=^ = 16'^ ■ 32 at /5 = 6 agrees only roughly with other estimates (within errors): the 
numbers are systematically higher than the central value and have larger statistical errors. It 
has been argued that these deviations can be partly related to the lattice Gribov ambiguity 
and that the breaking of translational invariance in the time direction may be responsible 
for an enhancement of systematic errors from finite volume effects. 



8 Gauge fixing and Confinement 

In the past decades, many explanations of the QCD confinement mechanism have been pro- 
posed, most of which share the feature that topological excitations of the vacuum play a ma- 
jor role. Depending on the underlying scenario, the excitations giving rise to confinement are 
thought to be magnetic monopoles, instantons, dyons, centre vortices, etc.. These pictures 



include, among others, the dual superconductor picture of confinement |7^, |Tl|, |7^, |7^, |7^, |75 
and the center vortex model 0, [f^, |T|, |7S|, 0. 

Different features of the infrared collective degrees of freedom dominating these two mod- 
els can be identified and isolated in different gauges: the Abelian gauges and the Center 
projection gauges respectively. These procedures end up with the SU (N) link variables pro- 
jected as close as possible to the elements of U{N) and Z{N) respectively. These procedures 
are explicitly gauge dependent and therefore it is relevant to analyze the projection-physics 



20 



dependence on the details of the gauge fixing procedure. A brief description of the more 
favored abehan and center projections, will be given in the following. 

Reviews on the confinement studies and discussions on the implications of the lattice calcu- 
lations for the question whether it is really monopoles or vortices that drive the confining 
physics or these idea are not necessarily exclusive can be found in Refs. |]8T|, |82|, p3 



8.1 Maximal Abelian Gauge 

In the scenario of 't Hooft and Mandelstam the QCD vacuum state behaves like a magnetic 
superconductor. A dual Meissner effect is believed to be responsible for the formation of thin 
string-like chromo-electric fiux tubes between quarks in SU{N) Yang-Mills theories. This 
confinement mecahnism has been established indeed in compact QED ||8^. The disorder 



of the related topological objects, magnetic monopoles, gives rise to an area law for large 
Wilson loops and, thus, leads to a confining potential. 

Nonperturbative investigations of this conjecture became possible after formulating the lat- 
tice version of 't Hooft's Maximal Abelian Gauge projection (MAG) [0. The idea is 
to partially fix gauge degrees of freedom such that the maximal abelian (Cartan) subgroup 
(?7^~^(1) for SU{N) gauge group) remains unbroken. Lattice simulations |^8|, |89[ have in- 
deed demonstrated MAG to be very suitable for investigations of SU (2) abelian projections. 
In the case of SU{2) gauge theory, fixing MAG on the lattice amounts to maximizing the 
functional 

Fu[G] = Y,Tr {a;U^{xf^^^a^U,{xYG{x)) (58) 

with respect to local gauge transformations G{x). Condition Eq. (pSf ) fixes G{x) only up to 
multiplications G{x) W{x)G{x) with W{x) = exp{ia{x)T3) , T3 = 0-3/2,— 2n < a{x) < 
Sii, i.e. G{x) G SU{2)/U{1). An SU{2) subgroup method ||90| can be used to perform the 
maximization of the diagonal components of the gauge fields with respect to the off-diagonals. 
The matrix diagonalization can be performed iteratively using local gauge transformations 
PT| and overrelaxation can be used in the gauge fixing procedure. 



After that configuration has been transformed to satisfy the MAG condition, the coset 
decomposition: 

f/^(x) = C^ix)V,,ix) (59) 

is performed, where V^{x) = exp{i^{x)T3), — Stt < <P{x) < 2-71, transforms like a (neutral) 
gauge field and C^(x) like a charged matter field with respect to transformations within the 
residual abehan subgroup 

V^{x) W{x)Vf,{x)W\x + fi). C^{x) W{x)C^{x)W\x + /i) (60) 

Quark fields are also charged with respect to such U{1) transformations. The abelian lattice 
gauge field V^{x) constitutes an abelian projected configuration. 

The SU{2) action of the original gauge theory can be decomposed into a f/(l) pure gauge 
action, a term describing interactions of the U (1) gauge fields with charged fields, i.e. the off- 
diagonal components, and a self-interaction term of those charged fields p2[- Maximizing 
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the diagonal components of all gauge fields with respect to the off-diagonal components 
amounts to enhancing the effect of the pure U{1) gauge part in comparison with those 
contributions containing interactions with charged fields. The MAG projection (and various 
abelian projections) might enhance the importance of the U{1) degrees of freedom and the 
V^(x) abelian gauge field can be used to investigate: Creutz ratios and Polyakov lines p3| , 

dual London relations 



monopoles densities [§8| |94 



74, 



expectation values of monopole 
creation operators |96[, disorder parameters relative to monopole condensation |]7^, etc.. 

Gribov ambiguities in the actual projection procedure on the lattice will be discussed in 
Section |Il|. In the continuum, the maximally Abelian gauge, its defining functional and the 
Gribov problem (the presence of Gribov copies is shown explicitly) in this gauge are reviewed 
and analyzed in depth in Ref. P7[]. 



8.2 Maximal Center Gauge 



The old idea about the role of the center vortices in confinement phenomena |Tg] has been 



revived recently with the use of lattice regularization. In particular, it has been argued that 



the center projection might provide a powerful tool to investigate this idea ||7^. The gauge 
dependent studies were done in center gauges leaving intact the center group local gauge 
invariance. It is believed that gauge dependent P-vortices (projected vortices) defined on the 
lattice plaquettes are able to locate thick gauge invariant center vortices and thus provide the 
essential evidence for the center vortex picture of confinement. After gauge fixing the link 
variables are projected onto the centre, i.e. they are replaced by the closest centre element. 
This procedure is in complete analogy with abelian projection in the abelian gauges and the 
center dominance (the analog of the abelian dominance) means that the projected string 
tension az{2) {<^u{i)) is very close to the nonabelian theory string tension <Jsu(2)- Maximal 
center gauges are defined in the lattice formulation of Yang-Mills theory by the requirement 
to choose link variables on the lattice as close to the center elements of the gauge group as 
the gauge freedom will allow. Attempts in this direction stemmed from the studies given in 
Ref. [^]. So far three different center gauges have been used in pratical computations: the 
direct maximal center gauge (described in the following), the indirect maximal center 
gauge |]78| and the laplacian center gauge ||98|, the simple center projection gauge fixing 



procedure has furtherly been proposed ||100| 



The direct maximal center gauge (DMC), widely used in SU{2) lattice confinement studies, 
is defined by the maximization of the functional: 

Fu[G] = {iTrU.ixf^A' = J2l {TraMxf^^^ + l) , (61) 



with respect to local gauge transformations G{x) ([l|). Condition (|6TD fixes the gauge up 
to Z{2) gauge transformation, and can be considered as the Landau gauge for adjoint 
representation. Any fixed configuration can be decomposed into Z{2) and coset parts: 
Ufj,{x) = Z^{x)V^{x), where Z^{x) = sign(Trf/^(x)). The plaquettes ^'^^^(x) constructed 
from the links Z^{x) have values ±1. The P-vortices (which form closed surfaces in 4D 
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space) are made from the plaquettes, dual to plaquettes with Z^^{x) = —1. The Center 
projection procedure is in complete analogy with Abelian projection in the MAG. 

The center gauges have been implemented generally for SU(2) |^ |79|, ^ |101| , |102| , 
|103| , |104| , |105| , |106|| . An alghorithm implementing DMC in SU{N) lattice gauge theory is 



proposed and checked on SU{3) vortex-like configurations ||107|| in Ref . ||108 



Nevertheless recent alarming results 12^, |109|| on the Gribov problem severity (see Section |TTD 
in the DMC procedure cast some doubts on the physical meaning of P-vortices in this gauge. 

The continuum analog of the maximum center gauge, in particular of the Polyakov gauge in 
which the Polyakov loop has diagonalized, can be found in Ref. | ]11U[ |. 



8.3 Issues on string tension 

A matter worthy of note is the measure of the string tension extracted from Wilson loops 
constructed from abelian and center projected link variables. The observation that the 
reduced theories reveal the full string tension (i.e the abelian or the center dominance) 
nurtures the conjectures that those degrees of freedom give rise to confinement. A series of 



studies ||111|| of numerical investigations in SU (2) lattice gauge theory has established that 
the abelian projection obtained with the MAG fixing indeed accounts for most of the string 
tension. Recent calculations of this quantity for SU{3) can be found in Refs. | 112 | | 113 | [ 114 
which results are consistent with the values quoted in the literature ||115| , |116| |. In particular 
in Ref. ||114|| a stochastic gauge fixing method which interpolates between the MAG and no 
gauge fixing is developed. The heavy quark potentials derived from Abelian, monopole and 
photon contributions is studied. For Abelian and monopole contribution it is observed that 
the confinement force is essentially independent of the gauge parameter. On the contrary the 
Gribov ambiguity (see Section |1T]) seems to influence severely the string tension in the case of 
the MCG. Recent studies ||109| , ^ contradict previous numerical simulations demonstrating 
that the entire asymptotic string tension was due to vortex-induced fluctuations of the Wilson 
loop [|T7 . 



9 Gauge Fixing in the Langevin Scheme 

In this Section we will discuss the use of gauge fixing in algorithms in which the field config- 
urations are generated using a discretized Langevin equation. This numerical technique has 
been adopted to implement the so-called difference method ||118| , |119| , p.20|| whose relevance, 
both for gauge and spin systems, has been widely acknowledged (see for example Ref. | ]ll(j| ]). 

In the Monte Carlo approach, correlation functions are obtained as expectation values of 
operators over a suitable number of uncorrelated configurations. As they are generally 
expressed as differences between similar numbers, they can be affected by large statistical 
fluctuations. The difference method is an interesting attempt to override this effect to get 
more accurate results. This method is based on the idea of perturbing the system far away 
from equilibrium in a limitated space-time region and measuring the decay of the correlations 
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from such zone by using specific properties of the dynamical updating algorithm. In this 
procedure one computes the differences between the matrix element values evaluated on 
a perturbed configuration and on the unperturbed ones. If one is able to keep the two 
configurations very close each other, a coherent cancellation of statistical errors between the 
two highly correlated stochastic processes follows in the difference. To describe the difference 
method let us consider two gauge systems K and K'{i) on two lattices of identical size M^xL; 
K is kept at inverse square coupling /3; K'{i) has a "time slice" i {i = 1,2,3...L) where the 
inverse square coupling takes the value [3 + 513. We define -E(j) as the average energy of 
the time slice j for the system K., and W{i,j) as the average energy of time slice j for the 



perturbed system K'{i). Then it is possible to show ||119|| that the connected correlation 



function at distance d of the energy operator C{d) ( using Wilson action) is 

for any i,j such that d = \j — i\. In order to avoid noisy correlations the two sets of 
configurations must be similar (we are forgetting for a moment about gauge invariance) and 
this happens only if the method used to generate the configurations is continuous in the f3 
variable. 



The Langevin I12TI , \12% [Tl8| , ITTgl , |l2l update scheme is a well known example of continuous 



algorithm. Here we will give some details of this scheme; it will be useful also in the following 
to define the numerical stochastic perturbation theory. If S is the action, the gauge field 
configurations can be obtained as solution of the following Langevin equation: 

UL{t) = ^^+VL{t), (63) 
where riL{t) is a gaussian noise with autocorrelation: 

{VL{t)VL'{t'))=6LL'S{t-t') . (64) 

Notice that the "time" t occurring in Eq. (^) has nothing to do with the physical euclidean 
time. In fact, it is the evolution time of the differential equation dynamics adopted to 
formulate Langevin algorithm in lattice gauge theory. After a certain lapse of time, the 
system will become representative of the Boltzmann distribution; therefore the Langevin 
equation provides a way to generate this distribution analytically and numerically. 

Some extra care is needed to discretize the Langevin equation as new parameters are needed; 
they have to be tuned in order to get good performances. Just to be concrete, we sketch 
here the implementation of the Langevin dynamics suitable for lattice simulations, from 
Ref . ||124|| . A single Langevin step is given by a sweep of the lattice where each link variable 



is updated according the rule 

U^{x) -> U'^ix) = e-^-(^) U^ix) (65) 



where FJx) is given by 



F^{X) = ^ Y^^Up - Ul)\traceless + VeH^{x) (66) 
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Here e is the Langevin time step; the sum over P means that gets contributions from all 
oriented plaquettes which include the link /i at x. Finally H^{x) is extracted from a standard 
(antihermitian, traceless) Gaussian matrices ensemble. 



9.1 The Chameleon Gauge 

The effectiveness of the difference method implemented by the Langevin dynamics depends 
on the evolution trajectories in the phase space of the two systems, unperturbed and per- 
turbed: they should be as close as possible. For a gauge model, the additional degrees of 
freedom create an additional complication. The gauge part performs a random walk in phase 
space, and tends to separate the two trajectories of the Langevin dynamics ||118| , |119| . A 



simple gauge fixing (for example putting all the time-like gauge variables = 1) turns out to 



slow down the dynamics and makes the method impractical ||119|| . Also the introduction in 



the Langevin equation of a magnetic field term ||119|| , which would ensure a smooth, partial 



gauge fixing, does not turn out to be successful. Trajectories in phase space diverge quite 
soon and there is a very unpleasant slow drift of the energy. 

To overcome these problems a peculiar gauge fixing called chameleon gauge was proposed ||125| 



The gauge is fixed in such a way that two fields are as similar as possible. This procedure 
does indeed keep the gauge part of two systems as close as possible reducing the rate of diver- 
gence of the two trajectories in phase space without introducing any sizable slowing down 
in the observable dynamics like the energy. After each full lattice sweep of the Langevin 
updates of link variables U (unperturbed fields) and V (perturbed fields), a gauge fixing is 
performed on the U configuration. One maximizes the quantity: 

Fu[G] = Y,U'',{x)Vi{x) (67) 

where x runs over the lattice sites and /i over the 4 directions, with respect to gauge transfor- 
mations G{x). The gauge fixing performed at each step has a sizable effect on the evolution 
of the system as the corresponding Eq. (^) is not gauge invariant. The quantity 

:f^Y,{^-U,{x)Vl{x)) (68) 

measures the gauge fixing quality, small values of T mean good gauge fixing. Correlations 
functions for the SU{2) 0^^ glueball mass measured on the chameleon gauged configurations 
are far less noisy (a factor of order 5) than with other methods. It is also quite interesting to 
note that the breakdown of the correlation functions is always signaled by a sudden growth 
of the quantity, i.e. by the collapse of the quality of the gauge fixing one is able to reach. In 
Fig. H it is shown the signal extracted at separation 1 by using the magnetic field method 
119|| , while in Fig. |^ it is shown the result obtained with chameleon gauge fixing. 
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Figure 6: The distance 1 correlation computed by using the magnetic field method. 
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Figure 7: The distance 1 correlation computed by using the chameleon gauge fixing. 
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9.2 Numerical Stochastic Perturbation Theory 



To compute large orders in the perturbative expansion of observables (therefore gauge invari- 
ant) in lattice field theories and also to the aim of overcoming some divergent fluctuations 
occurring in perturbative Langevin dynamics, the method of the numerical stochastic per- 
turbation theory (NSPT) has been proposed in Ref. | ]126| , |127|| . This method is implemented 



in the scheme of stochastic gauge fixing originally formulated in the continuum ||128|| and 
then extended to the lattice case ||129| . 



In this approach, perturbation theory is performed through a formal substitution of the ex- 
pansion {k is the perturbative order and g is the standard coupling in lattice gauge theories): 

U,{x)^J29'Ul''\x) (69) 

k 

in the Langevin equation (|65|) used for lattice simulations. The power expansion of the 
field induces a power expansion for every observable (y4^(x) included), then the perturbative 
expansions is usually computed as an average over the Langevin evolution. Even though the 
original motivation of the Langevin approach was to perform calculations in perturbation 
theory without fixing a gauge, it is known that some divergent fluctuations may plague high 
order terms (averaging to zero). A proposed way out is the technique of the stochastic gauge 
fixing. The underlying idea of this approach is the introduction of an attractive force in 
the Langevin equation in such a way that the field is attracted by the manifold defined by 
Landau gauge and that its norms are kept under control without affecting the observables. 
The implementation on the lattice consists in a gauge transformation, which is executed 
after each Langevin step, given by 

WL-.U^ix) e'"(")f/^(x)e-"'("+'=^) 

W{x) = a A^f,[U^{x) - Ul{x)]\\traceless 

A_^U,{x) = U,{x)-U,{x-e^) (70) 

One can prove | [L3CI| | that by doing this, the system gains a force that drives it towards the 
Landau gauge. By interleaving it to each Langevin step one obtains a sort of soft gauge 
fixing where Wl provides an additional drift which however does not modify the asymptotic 
probability distribution. After that, one has to expand in g the gauge fixing step Wl- this 
can be achieved with the same technique already developed for the unconstrained Langevin 
algorithm ||124|| . The value of the parameter a is chosen in such a way to minimize systematic 



errors. We report, for example, the term in g'^ of the plaquette measured in Fig. ^ without 
gauge fixing and in Fig. ^ with the above gauge fixing. The extension of NSPT with the 
adoption of the stochastic gauge fixing in a gauge non invariant context is possible with the 
caveat that also the gauge condition one wants to enforce has to be expanded as a series of 
conditions. To be definite, as the field is expanded as = A^/^^ + gA^f!''^ + g'^A^^^ + 
the form of Landau conditions one needs to impose is 

d,Ap = (71) 
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Figure 8: The fourth order in the expansion of the plaquette in a pure Langevin simulation. 




-1.5 

500 1000 1500 2000 

TIME 

Figure 9: The fourth order in the expansion of the plaquette in a simulation with stochastic 
gauge fixing. 
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for every order k (partial derivatives, as usual, are to be understood as finite difference 
operators). The implementation | |131|| goes as follows. The functional, as function of A^, to 
be maximized by the gauge transformation, is: 



X,fJ, 



(72) 



where an expansion for F4[iy] is induced by the expansion of the field A^. If the transfor- 
mation W is chosen as 



W{x) = ^g^w'^^\x) 

k>0 

^x;W(x) = , (0 < a < 1) 



(73) 
(74) 



then the extremum conditions for every N^'^^x) are recovered enforcing exactly eq. ([7l|). 
The strategy of the extension of NSPT to compute the expansion of gauge non invariant 



quantities |130, 131], is the following: 



i) let the system evolve (in the stochastic gauge fixing scheme) to get a thermalized 
configuration; 

ii) fix Landau gauge implementing the condition ( [7T| ) order by order and measure; 

iii) go back to i), i.e. let the system evolve until a decorrelated configuration is reached. 

The status of the method with respect to gauge fixed lattice QCD is revised and a first 
application to compact (scalar) QED is presented in Ref. [|130|| . A discussion about the 
convergence of the stochastic process towards the equilibrium, the expected fluctuations in 
the observables and the computations of quantities at a fixed (Landau) gauge in this frame 
can be found in Ref. ||132| | . A success of the NSPT extension has been the computation of the 
lattice SU{3) basic plaquette to order /3~^ ||127|| to actually verify the expected dominance of 



the leading infrared (IR) renormalon (a recent review on this item is in Ref.[ 133| ]) associated 
to a dimension four condensate. Recently the order has been computated ||134|| and 
then result is consistent both with the expected renormalon behaviour and with finite size 
effects on top of that. Another application of the NSPT method has been the computation 
of the perturbative expansion of the so called residual mass term in lattice heavy quark 
effective theory to order an|135|. 



10 Lattice Gauge Potential 

In this Section we will discuss the problems related to the ambiguities due to the lattice 
definition of the gauge potential A^. The root of the problem is that a unique, natural 
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definition of tlie potential on tlie lattice does not exist because in the Wilson discretization 
of gauge theories, the fundamental fields are the links [/^ which act as parallel transporters 
of the theory. Hence the lattice fields are derived quantities which tend to the continuum 
gluon field as the lattice spacing vanishes. As a consequence on the lattice it is possible 
to choose different definitions of formally equal up 0(a) terms and there is not any 
theoretical reason to prefer one or another. In quantum field theory this ambiguity is well 
understood because any pair of operators differing from each other by irrelevant terms, 
i.e. formally equal up to terms of order a, will tend to the same continuum operator, up 
to a constant, see for a general discussion Ref. ||136 |. In the Wilson's regularization the 



natural and most used definition of the 4-potential in terms of the links, ?7^, which represent 



the fundamental dynamical gluon variables is given by the standard relation ([2T|). This 
definition is certainly not the only possible and this ambiguity can create some problems in 
the discretizations of the continuum gauge fixing equations. Other definitions with analogue 
properties as, for instance 



A'Jx) = 



Traceless 



Aiag 



(75) 



which in fact differs from to the standard one (pi]) by terms of 0{a) that formally go to zero 
as a — > 0, have the same validity. 

Of course the requirement that the gluon fields Afj_{x) rotated in the Landau or Coulomb 
gauge satisfy the corresponding gauge conditions (see Section ^ is based on the gluon field 
definition, therefore different definitions generate different gauge fixings on the lattice 
which in the limit a ^ must correspond to the continuum gauge fixing. This feature, 
checked in perturbation theory, has been verified numerically at the non-perturbative level 
in Ref. | 137|| , where it has been shown that different definitions of the gluon field give rise to 



Green's functions proportional to each other, guaranteeing the uniqueness of the continuum 
gluon field. 

The relation between two A^ definitions can be expressed up to O(a^) terms in this way 



136 



A^{x)=C{go)A,{x). (76) 
Therefore for a Green's functions insertions the following ratio is expected to be a constant 

{...A'{x)...) 



{...A,{x)...) 



Cigo) ■ (77) 



This relation has been checked numerically on the lattice by measuring a set of Green's 
functions related to the gluon propagator for SU{3) in the Landau gauge with periodic 
boundary conditions computed with the insertion of the standard (^1]) and A'^ ([76|). In 
Fig. ([TOD {A[A'j) and the rescaled one Cf{go){AiAi) are shown, where 



(AA)(t) = -^5^5^Tr(A,(x,t)A,(y,0)) (7J 
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Figure 10: Comparison of the matrix elements of {A[A\){t) (crosses) and the rescaled 
(AiAi) • Cf{go) (open circles) as function of time for a set of 50 thermalized SU{3) configu- 
rations at (3 = 6.0 with a volume y • T = 8'^ • 16. The data have been slightly displaced in t for 
clarity, the errors are jacknife. From Ref. [137]. 



and the matrix elements {A[A[){t) are obtained replacing in the same form the alterna- 
tive definition (|7^). The remarkable agreement between these two quantities confirms the 
proportionality shown in Eq. {^^. 

From the algorithmical point of view, however, the various definitions are not interchange- 
able. In fact the parameter 6, monitoring the numerical behaviour of the gauge fixing 
algorithm (as explicated in Section as a function of lattice sweeps, is quite sensitive to 
different definitions. As shown in Fig. ([TI|) only the definition which appears in the 
F functional minimization, see Eq. (pOD, goes to zero. Note that the comparison reported 



in Fig. 11 is done on one configuration at time in order to check the gauge fixing quality 



while the comparison between the operators 6 and 9' must be done averaging them over the 
gauge fixed configurations of the thermalized set. Moreover the behavior shown in Fig. ([TT|) 
can be readily understood in the following way. The operator 6, defined in Eq. (P^D, and 
the operator {9 ), constructed with the same form with replaced by A'^, are computed 



in the lattice units taking the definition eq. (|2TD and Eq. ( [75D without the powers of a to 
the denominator. Then in the continuum variables 9=^1 d^^i'^^i^f^i^))'^ where V is the 
4- volume in physical units and analogously for 9' . Hence, while 9 vanishes configuration 
by configuration, as a consequence of the gauge fixing, 9 is proportional to {dfj^A'^)^, which 
has the vacuum quantum numbers and mixes with the identity. The expectation value of 
{d^A'^Y, therefore, diverges as ^ so that 9 will stay finite, as a ^ 0. 

As a matter of fact, the construction of lattice operators converging, as a ^ 0, to the funda- 
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Figure 11: Typical behaviour of 6 and O' vs gauge fixing sweeps at /? = 6.0 for a thermalized 



SU{'3>) configuration 8 • 16 from Ref. [137| 



mental continuum gauge fields, is affected, at the regularized lattice level, by an enormous 
redundancy due to the irrelevant terms. Although, if on the general field theoretical grounds 
the validity of such results is not unexpected, it is remarkable that it holds true also in this 
particular situation in which gauge fixing is naively performed, disregarding the problems 
related to the existence of lattice and continuum Gribov copies (see Section |ll]). 

The freedom to choose the lattice definition of can be used to build discretized functionals 



which lead to more efficient gauge-fixing algorithms. In Ref. ||138|| a new gauge fixing func- 
tional is proposed to remove the lattice discretization errors of order O(a^) to the Landau 
gauge condition. This improved scheme is used to fix the Landau gauge in the SU{2) lattice 
simulation of the gluon propagator in Ref. ||139|| . 

A Landau gauge fixing algorithm, using the exponential relation between link and gauge field 
is studied in Ref. 
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In Ref. [0 the gluon propagator is evaluated in different gauges and 
using different gluon field definitions on the lattice, corresponding to discretization errors of 
different orders. 



11 Lattice Gribov Copies 

In 1978 Gribov |^ discovered that for non-abelian gauge theories the usual linear gauge 
conditions does not fix in a unique way the gauge potential. In fact, it is possible to find 
different gauge potentials satisfying the gauge condition which are related each other by 
nontrivial gauge transformations. 

The presence of Gribov copies implies that the constraint of the naive gauge fixing is not 
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sufficient to remove all the degrees of freedom associated to the group of gauge transforma- 
tions. 

It is interesting to discuss the gauge fixing procedure in geometrical language, see for example 
Ref. ||141|| . The point where the gauge orbit, the curve described by field A'-^ as a function of 



the gauge transformation G, intercepts the plane defined by the gauge condition f{A'^) = 
represents a gauge-fixed potential. If the gauge orbit intersects in more than one point the 
plane f{A'^) = 0, each point represents a Gribov copy. It is possible to define a Hilbert 
norm of the gauge potential along the orbit as given in Eq. (^) . This definition has all the 
good properties of a norm and its values are able to distinguish among different local gauge 
transformations. Global gauge transformations does not change the norm and therefore they 
have to be considered in the same class of equivalence. The gauge transformation(s) capable 
to enforce the gauge condition can be found searching the stationary point (s) of the following 
functional (see also Eq. (|23|)): 

Fa{G) = -\\A^\\' = -J2 f d'xTrW,{x)f] . (79) 

In the case of minima of F the determinant of the Hessian matrix of the second derivatives 
of F is positive, this coincides with the Faddeev-Popov determinant (in the case of the 
Landau gauge, for example, the Faddeev-Popov operator FP is FP = —dD[A'^]). The 
set of the gauge potentials A'^ which are the minima of the corresponding F's defines the 
Gribov region Q which is known to be a convex region and its boundary dfl, where lowest 
eigenvalue of the Faddeev-Popov operator vanishes, is known as the (ffist) Gribov horizon. 
Note that the minima of F satisfy a more restricted gauge condition than the simple Landau 
gauge fixing which requires only the stationarity of F. This restriction, unfortunately, is 
not sufficient to solve the problem of Gribov copies because it is known that there can be 
multiple intersections also inside the Gribov region. The solution can be found by requiring a 
more restricted region of absolute minima, called the fundamental modular domain, contained 
inside the Gribov region. This scenario shows a possible attempt to solve the Gribov problem 
based on very interesting properties of the gauge potential topology. Unfortunately it is 
impossible to implement this constraint numerically. 

The studies of the Gribov problem on the lattice have demonstrated the existence of gauge 



fixing ambiguities both for abelian and non abelian theories for the first time in Refs. |p.42 



144] in the case of Coulomb and Landau gauges. It is remarkable that the scenario 



on the lattice seems to follow very closely the continuum one. In fact the lattice Gribov 
copies are found to correspond to different minima of the functional F so that they appear 
inside the lattice Gribov region. The solution could be to take the absolute minimum, which 
certainly exists in the compact theory, but this prescription is numerically hopeless. It 
must be noted also the corrispondence between the procedure and the formulas used on the 
lattice and in the continuum. Nevertheless, even if the phenomenological situation seems to 
be equivalent to the continuum one, it must be clearly understood that the Gribov copies 
are, at most or likely almost completely, determined by numerical lattice artefacts. On the 
other hand, in the continuum case the reason for the presence of Gribov copies is related to 
the deeper level of the theory, being usually connected with topological obstructions which 
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forbid the possibility to set up a smooth, diffential gauge condition vahd over all the lattice. 
Moreover, it is very difficult to make any connections between lattice Gribov copies and the 
continuum ones because the study of the continuum limit is operatively based on a chain of 
lattice simulations at different (3 values whose gauges change without any possibility to be 
controlled. Therefore the lattice studies of the Gribov ambiguity can have only a moderate 
influence on the theoretical uncertainties. Nevertheless the lattice Gribov copies may have 
a role in numerical simulations which therefore, must be kept under control. 

The generation of Gribov copies on the lattice is by now a standard method. Here we describe 
the procedure (mother and daughter method) which is schematically described in Fig. [12. 
Given a thermalized configuration M={?7}, where M is for mother, generated by a Monte 
Carlo simulation, one applies on it an ensemble of random, local gauge transformations. This 
is the cheap part of the procedure because the generation of random gauge transformations 
and the gauge rotations require a very short computer time. In such a way an ensemble 
of configurations gauge related to the mother are generated. Then, and this is quite 
expensive part in computer time, all these configurations are gauge fixed obtaining, at the 
end, the ensemble of the daughters dj. The final step consists of the analysis of the gauge- 
fixed ensemble, in order to understand whether all the members of the ensemble have been 
fixed to the same gauge configuration, or different minima (i.e. lattice Gribov copies) have 
appeared. In order to perform such a test, a good quantity to be measured is the final value 
of the functional F[G\ defined in Eq. (P^). In fact the F value is naturally gauge dependent 
and it is not affected by global gauge transformations which are to be considered gauge- 
equivalent, i.e. related to each other by global gauge transformations G{x) = G. Of course 
other gauge dependent forms |2^, |144| , |145|| may be adopted. 



To summarize: on the lattice the existence of many different minima of the functional F, 
not gauge equivalent, are called lattice Gribov copies and can be labelled with the value of 
the functional F itself. Of course it is unthinkable to succeed in reaching numerically the 
absolute minimum. The search of the F minima is at least as difficult as to find the lowest 
state of energy of a spin glass system with hamiltonian F. 

Besides the most known way of generating lattice Gribov copies based on the change of the 
gauge fixing process starting configurations by random gauge transformation, it is quite easy 
to produce Gribov copies on the lattice. As an example of that, we show in Fig. the lattice 
Gribov copies produced by the overrelaxation method when varying the overrelaxation 
parameter uo. In Fig. 0, different choices for the uo parameter, besides changing the the gauge 
fixing algorithm rate of convergence, lead to different Gribov copies even if starting from 
the same initial configuration. Of course there is no correlation between the convergence 
rate and the value of F associated with the particular lattice Gribov copy found. This can 
be qualitatively understood by recalling that non linear dynamical systems (and a gauge 
fixing algorithm is equivalent to such a system) often exhibit chaotic behaviour, in the sense 
that their evolution may depend dramatically both on the initial conditions (in this case 
the different random initial gauge rotations) and on the relevant parameter of the evolution 
equations (the overralaxation uj parameter). 

To remove the Gribov ambiguity it has been proposed to fix first to a uniquely defined axial 
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Figure 13: 6* as a function of the Landau gauge fixing sweeps for an SU{3) lattice for different 
values of the overrelaxation parameter u. The same letters 7,5) indicates curves with 



the same final F value, (from Ref. |36|) 
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But even if the starting gauge field is unique, one can still run into different 
local minima by using different gauge fixing algorithms. Moreover the axial gauge breaks 
rotational invariance, and some effect of this might be present in the local minimum that is 
favoured by this particular starting configuration. 



11.1 Gribov Copies and Measurements 



The numerical effects of lattice Gribov copies can be roughly divided into two categories: 



the distortion of a measurement and the lattice Gribov noise [17?]. The typical example of 
a distortion due to the existence of Gribov copies is the measure of the photon propagator 
in compact U{1) in the so called Coulomb phase. In this case the measure of the photon 
propagator as a function of the momentum, performed using the gauge fixing in the standard 
way, is affected by a not regular behaviour Jl^, [l^, jTof. This problem has been associated 
with the distortional effects due to the Gribov copies. In fact, after having chosen the gauge 
fixed configurations nearest to the minimum of the gauge functional, the photon propagator 
became a smooth, regular momentum function. More recent studies ||148|| show the details 
of the Gribov copies dynamics within the Lorentz gauge and provide a practical procedure 
to eliminate their effects in compact U{1) physics. In fact, in the Coulomb phase, it is 



shown that, apart from double Dirac sheets | 149 |, all Gribov copies originate mainly from 

The removal of these modes turns out 



m\ 



the zero-momemtum modes of the gauge field 
to be necessary for reaching the absolute maximum of the gauge functional. 

It is well known that Gribov copies have a relevant role in the lattice numerical studies on the 
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confinement vortex picture. A local iterative minimization in the direct maximal center gauge 
(discussed in the Section |8.2| ) selects every possible minima and each of these Gribov copies 
has its own set of P-vortices, which may show different properties: numerical demonstration 
of the Gribov copies effects has been shown in Ref. PT|. Starting from configurations fixed 



in the Lorentz gauge and then going to the maximal center gauge fixing one obtains, on 
average, a maximun higher than starting from random gauge and then applying the same 
gauge fixing procedure. The resulting picture in the two cases is dramatically different. In 
the latter case the very close numerical agreement of the center projected string tension az{2) 
of [^, |7S| with its currently accepted value 'unprojected' asu(2) is reproduced. In the former 
case there is essentially complete loss of information about the string tension value. 

Moreover, careful studies of the gauge copies problem for the direct center projection in 



SU{2) have been performed in Ref. ||109|| , indicating that this gauge is not appropriate for the 
investigation of the center vortices and that the Gribov copies problem is more severe than it 
was thought before. In ||109|| , following a procedure defined in Ref. ||151|| , a gauge dependent 



quantity X is computed on the gauge copy corresponding to the highest maximum of the 
functional F, Fmax{Ncop) after having generated Ncop gauge equivalent copies, for a given 
starting configuration. Averaging over statistically independent gauge field configurations 
and varying Ncop the function X{Ncop) is obtained and extrapolated to Ncop ^ oo limit. 
This procedue was built up in order to approach the global maximum as close as possible. 
The effect of the Gribov copies is clearly shown drawing the values of the Z(2)-projected 
Creutz ratios Xz{2){I) as a function of Ncop for (3 = 2.5, = 16^ plotted in Fig.|l^. Xz{2)il) 
is defined through the projected Wilson loops Wz{2)iC) = exp{i-n C{a p , C)} where C{<Jp, C) 
is the 4D linking number of the closed surface, ap, formed by P- vortex and closed loop C 

The Creutz ratios Xz{2){I) are used to estimate the projected string tension az(2)- An update 
of the " drama of Gribov copies" on the center vortices studies can be found in Ref. pOj . 



In this paper the disagreement between the projected string tension az{2) in the direct 
maximal center gauge and the physical string tension (Tsu{2) is demonstrated. The Laplacian 
center gauge ||^, appearing naturally as an extension of Laplacian Abelian gauge, is another 



attempt to get rid of Gribov ambiguity. 

It is interesting to note that in the measurements of the gluon propagator in SU{3) there 
is no signal about the Gribov ambiguity (for a recent review see Ref. 0). The numerical 
simulations are performed in the Landau gauge and the various authors claim that the effects 
of Gribov copies do not affect the measurement ||152|| . Anyway, in the more regular case in 
which there is no distortion due to Gribov copies, there should be an increase of the numerical 
fluctuations due to the incomplete gauge fixing associated with the copies. 

An attempt to study the properties of this noise has been done in Ref. ||153| , |154| , |64| taking as 



an example the measurement of the lattice axial current Z^- This quantity is particularly well 
suited to the Gribov fluctuations study. In fact Za is a gauge independent quantity which can 
be obtained from chiral Ward identities in two distinct methods: a gauge independent one, 
based on the matrix elements between hadronic states, and a gauge dependent one, based 
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Figure 14: The dependance of the Creutz ratios Xz(2)(-^) on the number of gauge copies N^op 
for /? = 2.5 as obtained in Ref. ||1U9| 



on the matrix elements between quark states in the Landau gauge. In the intermediate steps 
of the numerical computation, the second procedure takes into account gauge dependent 
matrix elements potentially subjected to the Gribov noise. Hence, there is an explicitly gauge 
invariant estimate of which is free of Gribov noise and which can be directly compared 
to the gauge dependent, Gribov affected, estimate. The results of the analysis ( ||153| , |154|| ) 
can be summarized in the following way: 

• there is a clear evidence of residual gauge freedom associated with lattice Gribov copies; 



the lattice Gribov noise is not separable from the statistical uncertainty of the Monte 
Carlo method. 



The global effect is not dramatic because the Za value obtained with the gauge dependent 
methods is close to the gauge independent evaluation and the jacknife errors are comparable. 

The influence of the Gribov copies has been studied also in the Coulomb gauge on the smeared 
correlation functions that are involved in the B meson decay constant determination. The 



B meson physics on the lattice involving gauge dependent operators like "smeared" 
sources for quark correlation functions, is a suitable sector for an investigation of the Gribov 
ambiguity. In fact, if the gauge condition implemented numerically (e.g. the Coulomb 
gauge) does not correspond to a complete gauge fixing, in the sense that the gauge fixing 
algorithm may converge randomly to any configuration in a set of Gribov copies , then the 
value of the operators will depend on which copy gets selected by the algorithm. Then such 
residual gauge freedom acts as a source of statistical noise in the Monte Carlo average of 
those physical quantities that are extracted from gauges dependent quantities. 
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It must be noted that the Coulomb gauge condition ([T8|) is implemented indipendently on 
each timeslice, given a starting link configuration. This is because the links in the time 
direction which connect adjacent timeslices do not appear in the definitions of In 
other words, each timeslice of a given configuration is endowed with its own scenario of 
Gribov copies, so that the pattern of their occurrence in the Coulom gauge is richer than 
in the Landau gauge, in which they are defined globally on the whole lattice. One finds 
that the residual gauge freedom associated to Gribov ambiguity induces observable noise 
effects, though at the level of numerical accurary of considered simulation these effects are 
not relevant to the final determination of Jb- The results in Ref. |TH], obtained from SU{3) 



lattice configurations generated on a 10^ x 20 lattice indicate that such effects may become 
important on bigger lattices. In fact, increasing the lattice size would typically reduce 
the standard statistical noise, allowing, in principle, to perform measurements from large 
time separations, but then the gauge noise effect (which is not expected to disappear in 
the continuum limit) may become a major source of fluctuations and may provide a relevant 
contribution to the error bars on the physical quantities evaluated from smeared correlations. 

Many numerical studies of Gribov ambiguities on abelian observables on the lattice have 
been performed in the last years. In the studies of the dual superconductor hypothesis of 
conflnement one mainly uses the lattice version of 't Hooft maximally Abelian gauge |57 



sketched in Section |8.1| . The effect of Gribov copies is clearly seen but the gauge flxing seems 
to be under control [ |151| , |155| , |15(j| , P7| . In particular, a careful study of the Gribov ambiguities 



in dual superconductor scenario has been performed in Ref. ||151|| . A new effective algorithm 



(simulated annealing ||155| | complemented with overrelaxation) flxing maximal abelian gauge 



to reduce gauge flxing ambiguities with respect to the standard overrelaxation algorithm and 
a numerical procedure to estimate the remaining gauge flxing uncertainties are used in this 
analysis. It is important to emphasize that, using a new method to assess uncertainties due to 
the incomplete gauge flxing, a procedure has been suggested to extrapolate values obtained 
on local maxima to the absolute maximum; the accuracy of all computations is limited by 



the statistical error on the biases. The investigation of Ref. ||151|| revealed that the effect 
of gauge copies cannot be neglected with respect to the statistical errors. Nevertheless, the 
proposed algorithm does reduce the variance of observables versus the gauge copies in a 
sizeable way and yields larger values of the functional to be maximized. Then the results 
for the string tension show that its abelian part accounts for 92 % of the conflnement part 
in the static potential (on a 32^ lattice at (3 = 2.5115). 



By analogy with the Laplacian gauge flxing a Laplacian Abelian gauge ||157|| (LAG), for the 



abelian projection, has been introduced. Also this gauge shares the smoothness properties 
with the Landau gauge but avoids lattice Gribov copies. This gauge flxing, which follows the 
construction of the Laplacian gauge, is based on the lowest-lying eigenvector of a covariant 
laplacian operator. A numerical study on abelian and monopole dominance in MAG and in 
LAG has been performed in Ref. [|158|| . An investigation and a comparison of the monopole 
structure of the SU{2) vacuum as seen in different gauges are also performed in Ref. Q . 
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Figure 15: Description of Ape's algorithm for smearing from Ref. ||168 

12 Gauge Dependent Smoothing Methods 

It is known that for the matrix elements of local operators the noise-to-signal ratio diverges 
as the continuum limit is approached (as discussed for example in Ref.[ [I59|| ). In particular, to 
check the scaling of physical quantitites, their dependence over /3 at the critical point of the 
theory has to be computed . Then it is necessary to go to higher (3 values, despite the fact that 
the lattice spacing shrinks when —* 0. Operators of a fixed lattice extent become smaller 
in physical size, and their projection over the ground state goes to zero. A possible trick in 
order to get a non-vanishing signal can be to resort to non local operators. Moreover, some 
kind of smoothing procedure is needed in order to study on the lattice the monopole-related 
observables and to determine the global topological charge in an unambiguous way ||160|| . The 



smoothing procedures kill part of the gauge field degrees of freedom in the process and these 
algorithms are not full gauge fixing methods. 

Several methods have been introduced to remove the unphysical short-distance fluctuations: 
cooling | p.61| , |162[ , smearing ||12CI| , |163| , |164|| , fuzzy loops | |165| , |166| , |167|| and other renormal- 



ization group smoothing methods (see for example Ref. [0 and reference therein). Here we 
will give just a brief description of smearing and fuzzy loops procedures. For the smearing 
we will discuss the gauge dependent procedure and for the cooling methods we send back to 
the quoted literature. 



The smearing procedure as originally proposed for SU (3) in ||120|| consists in the construction 
of correlation functions for operators which are a functional of the field smeared in space 
and not in time. Then smeared operators, usually used for the QCD spectroscopy and 
phenomenology on the lattice depend on the gauge and then it is necessary to fix the gauge 
before their measure. For each link of a configuration the product of the other three links 
defining a plaquette is considered , then these products are summed over the four choices of 
plaquettes orthogonal to the time axis; the resulting matrix, projected on the gauge group, 
is the new link variable. A graphic description ||168|| of the procedure is given in Fig.|15. 
The value of the smearing coefficient e is tuned in order to optimize the performance of 



the method [ 163 |. This procedure can be iterated to obtain a field more and more smeared 
in space at fixed time. Then one can have a set of operators, to implement a variational 
calculation, allowing to decide if a given euclidean time distance is asymptotic. Otherwise, 
as the amplitude of different exponential contributions depends on the operator, the ground 
state mass estimate can depend on the choice of the operator. In the momentum space the 
effect of the smearing is equivalent to the application of the factor exp(— A;^) (as it is shown 
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in Ref. |p,63|| for the example of a scalar field). This consideration makes more explicit the 
role of of short wavelength fluctuations suppressor taken by the smearing. Variations of the 
APE-style smearing have also been proposed in Refs. ||169| , |170 . 

The procedure of smearing may be used to improve any local operator involving link vari- 
ables. Smearing methods to improve lattice estimators have been already widely employed 
in the study of long distance correlations, such as large Wilson loops and hadron source 
operators. Just to give some explicit example we cite the improved lattice topological charge 
density operators, the hybrid (gluons as constituents) and glueball source operators. The 



first, constructed by a smearing-like procedure ||171| , |172| , |173|| , bears a better statistical 
behaviour as estimator of the topological density on the lattice both of SU{2) Yang-Mills 
theory and of the full QCD. For the second, it is necessary to use "smeared" or "fuzzed" 
links for the gluon parts of the operator to get good overlap with the ground state particles 
|174|| . To increase the overlap of hadron operators with their ground states the creation 
of a smeared quark source has been proposed in Ref. ||175| , p.76|| . This source, defined as 
S{x) = {D"^ + w^sc)m '^^>o where the parameter m^c was tuned to give the quark source an 
r.m.s. radius of about 3 lattice spacing, has been used for the the construction of smeared 
quark propagators. 

A smearing procedure, described in Ref. | |1 77| | , has been proposed to overcome the problem 
of Gribov copies on SU(2) lattices ||178|| . Following this procedure, for each link variable [/^ 
the sum i?^ of the 6 connecting staples is computed and is replaced with the combination: 



m = {U^ + wR^)/il + 6w) 



IV. 



where w is an adjustable parameter and is reunitarized. The procedure can be iterated 
many times; it turns out numerically that there is a critical value Wc of w below which 
the average plaquette goes to 1 (completely frozen lattice) for large values of the number 
of the smearing steps. The method is based on the following observation: when a large 
number of smearing steps are performed, the so called "trivial orbit" can be approached 
obtaining an unique Landau gauge configuration. Applying this gauge transformation to 
the original configuration, a unique starting point on the physical orbit is reached being the 
gauge path from this point to Landau point unambiguous. The smearing is stopped when 
a sufficiently frozen lattice is obtained. The method has been tested in SU{2) theory at 
/3 = 2.00,1.75,1.50. 

A different smoothing method, the fuzzy loops procedure, is depicted [168| in Fig. IB. In 
principle the coefficients are tunable even if ai = a2 is chosen ||167|| to normalize the 
fuzzy link back into SU{3). A consequence of fuzzy technique is that links paths grow 
exponentially during the iterations and there are 2^ times fewer fuzzy links than original links. 
Then, after the iteration, a simple loop of fuzzy links is a complicated linear combination 
of loops of original links. Finally elementary loops of fuzzy links are quite non local when 
expressed in terms of the original links. We note that the fuzzy procedure is inspired by the 



Monte Carlo renormalization group I^S], p.79| , |180|| methods involving factor-of-two blocking. 
Another renormalization group based smoothing method ||181|| is employed [0] on the lattice 
to investigate the monopole structure of the SU{2) vacuum. It is suitable to eliminate 
UV lattice artifacts from the monopole-related observables without destroing the confining 
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Figure 16: Description of Teper's algorithm for fuzzy links from Ref. [ |16^ . 



structure as other procedures finally do. 

Recently a new cooling approach ||182|| has been used as a gauge invariant low pass filter to 
extract physical information from noisy Monte Carlo configurations. In Ref. ||17CI|| Wilson 
action and topological charge are used to determine the relatives rates of standard coohng 
and smearing algorithms in pure SU{3) color gauge theory. 
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